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FAST QUANTUM MECHANICAL INITIAL STATE 
APPROXIMATION 

GOVERNMENT INTERESTS 

This application discloses an invention made with government support under 
Contract No. F30602-0 1-0523 awarded by US Air Force, Air Force Material 
Command Air Force Research Laboratory/IFKF. The government may have certain 
rights in the invention. 

FIELD OF THE INVENTION 

This invention relates to quantum computing and to methods and systems to 
efficiently calculate eigenvalues and eigenvectors Hennitian operators and quantum 
mechanical evolution operators with quantum computers and methods and systems to 
efficiently calculate an approximate quantum state to be used as an input in a quantum 
mechanical system. 

BACKGROUND 

Intuitively, quantum mechanical problems offer great potential for quantum 
computers to achieve large speedups over classical machines. An important problem 
of this kind is the approximation of an eigenvalue of a quantum mechanical operator. 
In a recent paper published in 1999 in Physical Review Letters (Vol 83, p. 5162) and 
hereby incorporated by reference, Abrams and Lloyd present a quantum method for 
doing this. Their method is exponentially faster than the best classical method, but 
requires a good approximation of the corresponding eigenvector as an input. 

There is currently a continuing need for a method and system for efficiently 
computing a good approximation of the eigenvector as an input to the Abrams and 
Lloyd quantum method. 
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There is also a continuing need for a method and system for efficiently 
computing a good approximation of a quantum state (not limited to eigenvectors) as 
an input to a quantum mechanical computer or computation. For example, one would 
like to compute an approximate input to the quantum simulation algorithm. The 
quantum simulation algorithm is described in the book Quantum Computation and 
Quantum Information, by M. A. Nielsen and I. L. Chuang, Cambridge University 
Press, Cambridge UK (2000). 

SUMMARY OF THE INVENTION 

The present invention is a system and method for use on a quantum computer 
to efficiently prepare the initial quantum state required by Abrams and Lloyd's 
eigenvalue approximation method. The system and method of the present invention is 
used to prepare a quantum register with an approximation of the eigenvector that is 
guaranteed to be sufficiently good to be used as input to the Abrams and Lloyd 
method. The present invention can be used when solving continuous Hermitian 
eigenproblems, e.g. the Schrodinger equation, on a discrete grid. 

Beginning with an eigenvector for a problem discretized on a coarse grid, the 
system of the present invention efficiently constructs, quantum mechanically, an 
approximation of the same eigenvector on a finer grid. This eigenvector 
approximation is suitable as the initial state for the eigenvalue estimation method of 
Abrams and Lloyd. 

Similarly beginning with a vector (i.e., a quantum state) for a continuous 
problem discretized on a coarse grid, the system of the present invention efficiently 
constructs, quantum mechanically, a vector (i.e., a state), which is an approximation 
to the corresponding vector on a finer grid. Our system efficiently extends a vector of 
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low dimension to one of high dimension, which is then presented as input to some 
quantum computation method, e.g., the quantum simulation algorithm. 

The features and advantages of the present invention will be more readily 
apparent and understood from the following detailed description of the invention, 
which should be understood in conjunction with the accompanying drawings 
appended to the end of the detailed description. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a chart illustrating the steps performed on the various quantum 
registers according to an embodiment of the present invention. 

Figure 2 is a chart illustrating the steps performed on the various quantum 
registers according to another embodiment of the invention. 

DETAILED DESCRIPTION 

For purposes of illustration only, and not to limit the scope of the present 
invention, the invention will be explained with reference to the embodiments of the 
invention indicated in the drawings. One skilled in the art would understand that the 
present invention is not limited to the specific examples disclosed and can be more 
generally applied to other initial state preparation methods and systems than those 
disclosed. 

The key component in the Abrams and Lloyd method is quantum phase 
estimation, which is a method for approximating an eigenvalue of a unitary matrix. 
Quantum phase estimation is also described in the above referenced book of Nielsen 
and Chuang. We give a brief outline of this method below. 

Let Q denote a 2"' x 2'" unitary matrix. We want to approximate a specific 

eigenvalue of Q. Phase estimation does this using the corresponding eigenvector as 

input. The Abrams and Lloyd method deals with the case when this eigenvector is not 
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known exactly. Referring to Figure 1, consider a quantum computer consisting of 
three registers 140, 150, and 160 with a total of b + m-^w qubits. The first b qubits in 
register 150 are all initially in the state |0>. The second register 140 with m qubits is 
initialized to some state | y/}, which must approximate the eigenvector in question 
sufficiently well, as will be seen. The last w qubits in register 160 are work qubits for 
temporary storage. The w qubits are not important in our discussion here, and we 
generally omit discussion of them below. 

Since Q is unitary and therefore normal, the state | y/} can be expanded with 
respect to eigenvectors of Q. Omitting discussion of the work qubits in register 160, 
the initial state of the algorithm is 

|0)|^) = 10)S*^»' <i) 

where \u} are the eigenvectors of Q. Placing the first register 150 in an equal 
superposition, using b Hadamard gates in step 170, transforms this state into 

1 ^i..^,^r...t (2) 



V2^ ^ V 



Next, powers of Q are applied in step 170 to create the state 




Since Q is imitary, its eigenvalues can be written as e^^ ^ ^ where (p „ ^ R- We 

can assume that (pu e [0,1) and consider the approximation of one of these phases 

instead of the approximation of one of the eigenvalues. Equation (3) is equal to 
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It is easily seen that the inverse Fourier transform performed in step 170 on the first 
register 150 creates the state 



i / 



where 



In step 180, a measurement of the first register 150 produces outcome j 190 with 
probability 

and the second register 140 will coUapse to the state 

represented by the register 200; register 210 contains the work qubits after the 
measurement 180 as known to one skilled in the art. 

We remark that for special case when the eigenvalues (p u can be represented 
exactly with fc-bits (i.e., 2^(pu is an integer), equation (5) simplifies to 
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When the eigenvalues are of tihis form and are distinct, a measurement in step 180 of 
the first register 150 will cause the second register 140 to collapse exactly onto the 
corresponding eigenvector in register 200. 

Recall that the system and method of the present invention are to achieve an 
approximation of the phase that corresponds to an eigenvector \u'} using a quantum 
computer, that the state \y/) is an approximation of this eigenvector, and that the 
eigenvalue is obtained from the value of the outcome j 190 by e^"" ' J^^""*" is of the form 
and approximates e^"^ For instance, one is often interested in the eigenvalue 
corresponding to the ground state or in low order eigenvalues. We define A((po, cpi) = 
min^G^ |.^ + ^1 - <Po\ }' Po.^i ^ ^ fractional part of the distance between cpo 

and <pi). Then a measurement of the first register produces an outcome from the set G 
= {j:A(j/2^,(p^.)<k/2\k>l } with probability 

and when A: = 1 the probability to obtain j such that a(j I2^,(p^)< is bounded 
from below by • I must be chosen in a way that this probability is large or 

preferably greater than 1/2, which implies that \dj\ has to be sufficiently large. For 
one embodiment of the present invention to obtain an approximation of (pu' with 
accuracy and probability at least \dj^ (l - s), equation (10) shows that the number 
of qubits h in the first register 150 must be 

6 



wo 2005/060400 



PCT/US2004/003074 



&=^+[^log(l + ~^^. (11) 

Quantxim phase estimation can be used as an efficient subroutine to find 
eigenvalues. Consider a Hennitian operator /f. The operator Gf^j = ^"'^^ is unitary 
and has the same eigenvectors as H. We assume that G can be implemented 
efficiently and, therefore, can be used as the unitary operator in the phase estimation 
algorithm. For example, when H is local, i.e., it can be written in the form , 

where each Hj acts only on a small number of qubits, then G can be implemented 
efficiently. However, locality is not a necessary condition for efficient 
implementation. Indeed, G can be efficiently implemented for a many-particle 
quantum mechanical system with a non-local H, One skilled in the art will 
understand that it is possible to implement G for a wide class of non-local 
Hamiltonians. 

The Hermitian eigenproblem described above is solved on a discrete grid. 
One embodiment of the present invention addresses the case in which the grid is 
extremely fine. Clearly, a fine grid requires a large vector for the representation of 
the initial state of the algorithm. In general, it may not be possible to efficiently 
prepare an arbitrary quantum state in a space with a large number of qubits. 
However, the present invention includes a method for the efficient preparation of an 
initial state. 

In one embodiment of the invention, the operator possesses an eigenvector for 
a coarse grid discretization of the problem, which was most likely obtained classically 
since the size of the problem is small, although one skilled in the art will understand 
an eigenvector obtained by any coarse method can be employed without diverging 
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from the scope of the invention. Using this eigenvector, we efficiently construct an 
approximation to the corresponding eigenvector for a fine grid discretization of the 
problem. We use this approximation as the initial state of the eigenvalue 
approximation algorithm. We describe our method for a one-dimensional continuous 
problem on the interval [0,1]. 

Let be a positive Hermitian operator, defined on a Hilbert space of smooth 
functions on [0,1]. Let vKO» = 1,2,..., denote the eigenfunctions of H, ordered 
according to the magnitude of the corresponding eigenvalues; and without loss of 
generality we assume that 

j\vu{^)f^^l^ (12) 

Suppose that is a discretization of H with grid size Hn = 1/(1 + AT). Let 
I [/^/^) ),A: = 0,l,...,iV-l,denotethenormalizedeigenvectors of ordered 

according to the magnitude of the corresponding eigenvalues. The expansion of the k- 
th eigenvector in the computational basis can be written as 

Let I y f ^ ) = 2 ((j + 1);^^ )\j ) be the sampled version of Vfc(-) at the 
discretization points. Consider problems such that the eigenvector of interest satisfies 
iK|L=supo^|v,(x)| = 0(l) and 
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where ^ > 0 is the order of convergence and j |z )|p = ^ 7=o l*''^/ f 

I X) = ^^.^Q ^ ) - For example, these conditions are satisfied when, for 
example, we are dealing with second order elliptic operators. 

Now, assume that the eigenvector |i7;[^''^^ of H^^ has been obtained 
classically. This vector is placed in a log iV^qubit register 110 (see Figure 1). ForiV= 
2 Wo, we construct an approximation |c/j[-^^^ of by appending s qubits in 

register 120, each qubit in the state |0) , to |i7j[^°^^ and then performing in step 130 a 
Hadamard transformation on each one of these s qubits in register 120, i.e. 



ps; (15) 



where f(j) = [j/^" J. The effect of /is to replicate the coordinates of \ul^'^)2' 

times. According to the present invention, [c/^^^^ is used as input to the eigenvalue 
and eigenvector approximation method. When the result of the method is measured 
1^1^^) will collapse onto a superposition of eigenvectors according to equation (8). 

The magnitude of the coefficient of jc/j^^^) in this superposition can be made 
arbitrarily close to one by appropriately choosing No. 
Consider two different expansions of U^^^"^ : 
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(16) 
(17) 



The first expansion is in the computational basis and the second is with respect to the 



eigenvectors //at. We call 



rewritten as 



d^^f} the probability of success. Equation (17) can be 



l9£jb (18) 

Taking norms on both sides and using (13) and (16) gives the inequality 



- 1-|<W' (IS) 
We will now bound (19) jfrom above, and thus the probability of failure. The 
definition of ^U^^^'j implies 



where EjJoKif = ^fe^ Zj^^'K^o/ = ro(h^J by (14). Applying fee 
triangle inequality, we get 
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The definition of \v^^^) and the fact that jv^j^ = 0{l) imply that 
IVj^") )||= ViV (l + 0{hif )) . Hence, the sum above is equal to 

Since Vft(-) is continuous with a bounded first derivative, we have that 

where J!^^^ = (7 + i)hf, and jCj.; = (/ U) + l)^Af„ J = 0,...,N-1. Clearly 

|-^2,i ~ I = ^{^^N„ ) • Using (22), (23) and the triangle inequality, we obtain fixim (21) 

that 

Hence, the probability of failure is bounded from above by o(a^o "^^^'^^^j • ^ 
depends only on the order of convergence to the continuous problem and the number 
of points in the classically solved small problem. We can select an No such that the 
probability of failure is less than Vi, no matter how much larger N is. By choosing a 
large N, we can make the discretization error arbitrarily small. Equation (24) implies 
that the probability of obtaining the eigenvalue with accuracy 2"^ is at least 

We remark that any classical numerical algorithm that computes an 

eigenvalue, satisfying a specific (nontrivial) property, of a X iV unitary matrix takes 

time f2(iV). For example, one may want to find the eigenvalue that corresponds to the 

ground state. This is true even if a matrix is sparse and regardless of whether the 
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algorithm is deterministic or randomized. It is merely a consequence of the fact that 
the algorithm needs to consider all the (nonzero) elements of the matrix, and there are 
at least Ci(N) such elements. Alternatively, in the restricted case when the matrix is 
diagonal finding one of its elements is a problem at least as hard as searching an 
unordered list. The lower bound for searching yields the lower bound in our case. 

In conclusion, our method provides a highly efficient preparation of initial 
states for eigenvalue approximation, requiring only a small number of Hadamard 
gates. Thus the method of Abrams and Lloyd, using the initial state prepared by the 
system and method of the present invention, computes the eigenvalue exponentially 
faster than any classical algorithm. The method of the invention can be generalized to 
higher dimensional continuous problems. 

In another embodiment of the invention, if we possess a vector that 
corresponds to a coarse discetization of a continuous problem then, under suitable 
conditions, we can efficiently extend it to a vector that approximates the 
corresponding vector (i.e., a state) of a fine discretization. Referring to Figure 2, we 
first place the original or given vector in register 310. Assuming that the vector has 
dimension No this register has log Nq qubits. For a AT = 2^No, we append to register 
310 s qubits, in the state | O) , in register 320. Then in step 330 we apply the 

Hadamard transform to the appended qubits. See Equation (15) and the explanation of 
the effect of the replicating function/. In register 340 we have the combination of the 
two registers 310, 320, register 340 containing the approximation corresponding to a 
vector (i.e., state) of dimension = 2^No. This requires log = log No +s qubits for its 
quantum mechanical representation. Step 350 represents a quantum mechanical 
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system using the approximation obtained in register 340. Step 360 represents the final 
state of the system 350. 

Having described the embodiments of the invention, it should be apparent that 
various combinations of embodiments may be made or modifications added thereto as 
is known to those skilled in the art without departing from the spirit and scope of the 
invention. 
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